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Abstract 

A systematic study of isotope chains in the rare-earth region is presented. For the chains 
^^^~^^^Nd, ^^^-^^°Sm, ^^^-^llGd, and ^^^'^HDy, energy levels, E2 transition rates, and two- 
neutron separation energies are described by using the most general (up to two-body terms) IBM 
Hamiltonian. For each isotope chain a general fit is performed in such a way that all parameters 
but one are kept fixed to describe the whole chain. In this region nuclei evolve from spherical to 
deformed shapes and a method based on catastrophe theory, in combination with a coherent state 
analysis to generate the IBM energy surfaces, is used to identify critical phase transition points. 

PACS numbers: 21.60.-n, 21.60.Fw 

Keywords: quantum phase transitions, catastrophe theory, interacting boson model 



I. INTRODUCTION 



Recently, a renewed interest in the study of quantum phase transitions in atomic nuclei 
has emerged A new class of symmetries that applies to systems localized 

at the critical points has been proposed. In particular, the "critical symmetry" E{5) 
has been suggested to describe critical points in the phase transition from spherical to 7- 
unstable shapes while X{5) 6] is designed to describe systems lying at the critical point 
in the transition from spherical to axially deformed systems. These are based originally on 
particular solutions of the Bohr-Mottelson differential equations, but are usually applied in 
the context of the Interacting Boson Model (IBM) Q , since the latter provides a simple but 
detailed framework in which first and second order phase transitions can be studied. In the 
IBM language, the symmetry E{5) corresponds to the critical point between the U{5) and 
0(6) symmetry limits, while the X{5) symmetry should describe the phase transition region 
between the U (5) and the SU (3) dynamical symmetries, although the connection is not a 
rigorous one. Very recently the 0(6) limit itself has also been proposed to correspond to a 
critical point 

Usually, the IBM analyses of phase transitions have been carried out using schematic 
Hamiltonians in which the transition from one phase to the other is governed by a single 
parameter. It is thus necessary to see how much these predictions vary when a more general 
Hamiltonian is used. Th e g lobal approach was first used by Castanos et al for the study 
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of series of isotopes 
consistent Q formalism (CQF) 



. An alternative procedure is provided by the use of the 
13l |. In this case, although the Hamiltonian is simpler than 



the general one, the main ingredients are included. Within this scheme a whole isotope 
chain is described in terms of few parameters that change smoothly from one isotope to 
the next. Because of the possible non-uniqueness of such nucleus by nucleus fits and the 
restricted parameter space, it is important to study under what circumstances the prediction 
of the location of critical points in a phase transition is robust. In this paper we follow 



Refs. 
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I4L I5I and use a more general one- and two-body IBM Hamiltonian to 



obtain the model parameters from a fit to energy levels of chains of isotopes. In this way a set 
of fixed parameters, with the exception of one that varies from isotope to isotope, is obtained 
for each isotope chain and the transition phase can be studied in the general model space. 
The fit to a large data set in many nuclei diminishes the uncertainties in the parameter 
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determination. A possible problem arising from working with such a general Hamiltonian, 
however, is the difficulty in determining the position of the critical points. Fortunately, the 



methods of catastrophe theory [3] allow the definition of the essential parameters needed 
to classify the shape and stability of the energy surface jl^ . isl ]. 

In this paper we analyze diverse spectroscopic properties of several isotope chains in the 
rare-earth region, in which shape transition from spherical to deformed shapes is observed. 
We combine this study with a coherent-state analysis and with catastrophe theory in order 
to localize the critical points and test the X(5) predictions. Since the introduction of the 
E{5) and X{5) symmetries, only a small number of candidates 
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3, HE, 
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2d, 3,0 



have been proposed as possible realizations of such critical point symmetries. In this paper 



we show t 
approach 



lat the critical points can be clearly identified by means of a general theoretical 



1$. 



The paper is structured as follows. In section |n] we present the IBM Hamiltonian used. 
In section IIIII the results of the fits made for the different isotope chains are presented. 
Comparisons of the theoretical results with the experimental data for excitation energies, 
E2 transition rates and two- neutron separation energies are shown. In section HVl the intrinsic 
state formalism is used to generate the energy surfaces produced by the parameters obtained 
in the preceding section. In addition, the location of the critical point in the shape transition 
for each isotope chain is identified by using catastrophe theory. Also in this section, the 
alternative description provided by the CQF for the rare-earth region is briefly discussed. 
Finally, section El is devoted to summarize and to present our conclusions. 



II. IBM DESCRIPTION 



In this work we use the interacting boson model (IBM) to study in a systematic way the 
properties of the low-lying nuclear collective states in several even-even isotope chains in 
the rare-earth region. The building blocks of the model are bosons with angular momentum 
L = {s bosons) and L = 2 {d bosons). The dynamical algebra of the model is U{6). 
Therefore, every dynamical operator, such as the Hamiltonian or the transition operators, 
can be written in terms of the generators of the latter algebra. Usually some restrictions 
are imposed on these operators, e.g. the Hamiltonian should be number conserving and 
rotational invariant, and in most cases it only includes up to two-body terms. 
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The most general (including up to two-body terms) IBM Hamiltonian, using the multi- 
polar form, can be written as 

H = AN + B ^^^~ + Sana + kqP^P 

+ KiL ■ L + K2Q ■ Q + K3T3 ■ fa + 1^4X4 ■ (1) 

where iV, and fid are the total boson number operator, and the d boson number operator, 
respectively and 

= \ {d^ -d^ - ■ s^), (2) 



T0(rf^ X df^\ (3) 

Q = (s^ X d + S X - ^(dt X J)(2), (4) 

fs = {S X df^\ (5) 

= (rft X df^\ (6) 

The symbol ■ stands for the scalar product, defined as T^-Tl = J2m{~^)^^ T^mTl-m where 
Tlm corresponds to the M component of the operator f^. The operator = (— l)™'7£-m 
(where 7 refers to s and d bosons) is introduced to ensure the correct tensorial character 
under spatial rotations. 

The first two terms in the Hamiltonian do not affect the spectra but only the binding 
energy. Therefore they can be removed from the Hamiltonian if only the excitation spectrum 
of the system is of interest. However, a complete description of both excitation and binding 
energies requires the use of the full Hamiltonian (Q). 

The electromagnetic transitions can also be analyzed in the framework of the IBM. In 
particular, in this work we will focus on E2 transitions. The most general E2 transition 
operator including up to one body terms can be written as 



- p „ 



{s^ X d + d^ X s)^^j + X{d^ X d)fj , (7) 



where Cg// is the boson effective charge and x is a structure parameter. 

Two-neutron separation energies (S2n) are also studied in the present work. This observ- 
able is defined as the difference in binding energy between an even-even isotope and the 
preceding even-even one, 

S2n = BE{N)-BE{N-l), (8) 



where corresponds to the total number of valence bosons. Note that if only the first two 
terms in (H)) are considered and A and B are assumed to be constant along the isotope chain, 
would be given by 

S2n = -{^-lB)-BN = A + BN. (9) 
For a detailed study of this property we refer to Ref. 



III. FITS 

In this section we analyze several isotope chains belonging to the rare-earth region using 
the most general IBM Hamiltonian, Eq. and E2 transition operator, Eq. (|7j). As an 
ansazt for each chain of isotopes we will assume a single Hamiltonian, and a single E2 
transition operator. All parameters in these operators are kept fixed for a given isotope 
chain, except for the single particle energy which is allowed to vary slightly from isotope 
to isotope. The way of fixing the best set of parameters in the Hamiltonian is to carry 
out a least-square fit procedure of the excitation energies of selected states (2]'', 4^, Q'^, 
^i, 0^, 23 , 43 , 2^, 3]*", and 4^) and the two neutron separation energies of all isotopes 
in each isotopic chain. Once the parameters in the Hamiltonian are obtained, the B{E2) 
transition probabilities 2^ Of, Af 2^, 2^ Of, 2^ Of, 0^ 2^, and 0^ 2t 
of the set of isotopes are used to fix Ce// and x by carrying out a least-square fit. The 
experimental data for excitation and binding e nerg ies and B(E2)'s have been taken from 



Refs. 
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Finally, it is worth noting that in 
Ref. |25| the Hamiltonian parameters were fixed just using the data for excitation energies 
and then A and B were adjusted to reproduce the experimental values of S2n- In this paper, 
since we are particularly interested in accurately describing the spectroscopic data associated 
to shape transitions, both, excitation and binding energies, are treated on an equal footing 
describing the shape transition, to determine the set of Hamiltonian parameters in Eq. (P). 

Tables in and ITTl summarize the parameters obtained for the Hamiltonian and E2 transition 
operator for each isotope chain. 

In figures ^ |21 El and |3 the systematics of experimental and calculated energies for the 
states included in the least-square procedure are presented in order to show the goodness 
of the fitting procedure. In figures El IHl 13 and |H1 the systematics of the experimental and 
calculated B{E2) values are compared. Finally, in figure IHl the experimental and calculated 
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values are shown. This is a fundamental magnitude for identifying a phase transition 
since it is directly related to the derivative of the energy surface. First order phase transitions 
are related with the appearance of a kink in the S2ri values. As shown in Fig. El the 
calculation matches the experimentally observed behavior. 

The analysis of the preceding figures for different observables and for several isotope 
chains shows that the present procedure is appropriate for systematic studies and confirms 
that it provides a simple framework to describe long chains of isotopes and detect possible 
phase transitions. 

An alternative approach to describe long chains of rare-earth nuclei is to use the CQF. 
The CQF Hamiltonian is 

H = ehd + kQ' ■ Q' , (10) 

with 

Q' = {s^ xd + d^ X sf^ + x{d^ X df^\ (11) 

For each nucleus the parameters e, n and x ^I's determined in order to fit the excitation 
energies and B(E2)'s. In particular in Ref. jQ the parameters of the Hamiltonian are 
calculated within the CQF framework with the ansazt that the strength of the quadrupole 
term of the Hamiltonian remains constant along a wide region of the mass table. As in the 
present paper they compare experimental data and theoretical values for excitation energies 
and B{E2) transition rates. Both methods provide a consistent description of the rare-earth 
region with a similar number of parameters as can be observed in Fig. ^1 and in table IIIII 
where the case of ^^^Sm is analyzed. Note that in the present work the results come from 
a global analysis, therefore the B{E2) transition rates are not normalized to the transition 
B[E2 : 2i Oi) in a particular isotope. If in table ITTTl the results are normalized so as to 
reproduce the observed value for B{E2 : 2 Oi") in ^^^Sm the results of this work and 
CQF are basically the same. 



IV. ENERGY SURFACES AND PHASE TRANSITIONS 

The studv of phase transitions in the IBM requires the use of the so called intrinsic-state 
formalism 0, ^ l^l, although other approaches can be used 0, |4^. This formalism is 
very useful to discuss phase transitions in finite systems because it provides a description of 
the behavior of a macroscopic system up to effects. To define the intrinsic, or coherent. 



state it is assumed that the dynamical behavior of the system can be described in terms of 
independent bosons ( "dressed bosons" ) moving in an average field . The ground state of 
the system is a condensate, |c), of bosons occupying the lowest-energy phonon state, Fj, 

\c) = ;^(rl)'^|o), (12) 



where 



^yYZ^ y ' /^cos7 4 + ^/5sin7(4 + c?-2)j (13) 
and (3 and 7 are variational parameters related with the shape variables in the geometrical 
collective model. The expectation value of the Hamiltonian in the intrinsic state p2|) pro- 
vides the energy surface of the system, E{N,f3,'y) = {c\H\c). The energy surface in terms 
of the parameters of the Hamiltonian ffl and the shape variables can be readily obtained 

0, 

NP"^ ( 9 7 9 



iV(iV-l) 

(1+/52)' 



J + /5'(- Y + 4 K^) + 2 72/^3 cos(3 7) 



where the terms which do not depend on /5 and/or 7 (corresponding to A and B in Eq. (|T|)) 
have not been included. 

The equilibrium values of the variational parameters /? and 7 are obtained by minimiza- 
tion of the ground state energy (c|H^|c). As mentioned above these parameters are related 
to the parameters of the Geometrical Collective Model and provide an image of the nuclear 
shape for a given IBM Hamiltonian. A spherical nucleus has a minimum in the energy 
surface at /3 = 0, while for a deformed one the energy surface has a minimum at a finite 
value of /? and 7 = (prolate nucleus) or 7 = 7r/3 (oblate nucleus). Finally, a 7-unstable 
nucleus corresponds to the case in which the energy surface has a minimum at a particular 
value of /? and is independent of the value of 7. The equilibrium values of /? and 7 are the 
order parameters to study the phase transition of the system, although in the case under 
consideration (IBM-1) only /5 has to be taken into account, since the minima in 7 are well 
defined. 

In Fig. ^2 the energy surfaces for the isotopes of the different isotope chains studied in 
this paper are plotted as a function of /5. The figure on the right is a zoom of the region 
close to /3 = 0. 
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The classification of pfiase transitions tliat we follow in this paper and that is followed 
traditionally in the IBM is the Ehrenfest classification In this context, the origin of a 
phase transition resides in the way the energy surface (their minima positions) is changing 
as a function of the control parameter that, in this work, is a combination of parameters 
of the Hamiltonian (see Eq. ((21]) )• First order phase transitions appear when there exists a 
discontinuity in the first derivative of the energy with respect to the control parameter. This 
discontinuity appears when two degenerate minima exist in the energy surface for two values 
of the order parameter (3. Second order phase transitions appear when the second derivative 
of the energy with respect to the control parameter displays a discontinuity. This happens 
when the energy surface presents a single minimum for /3 = and the surface satisfies the 
condition = 0. 

With the introduction of the E{5) and X{5) symmetries to describe phase transitional 
behavior, diverse attempts to identify nuclei that could be located at the critical points have 
been made. The theoretical approaches have been mainly performed with restricted IBM 
Hamiltonians. In particular, within the CQF, or other restricted Hamiltonians, the location 
of the critical point is obtained by imposing ^| = at /? = 0, where E is the energy surface 
P]. This condition leads to a fiat surface in a region of small values of (3, with a single 
minimum in the limit x = and two almost degenerate minima (one of them in /3 = 0) in 
the other cases. In the CQF approximation it can be said that (^l)^ ~ corresponds 
approximately to a "very flat energy surface" as happens for the E{5) and X{5) critical 
point models. Following this approach both ^^''Nd and ^^^Sm have been found to be close to 
critical. However, when studying a transitional region in which the lighter nuclei are spherical 
and the heavier are well deformed, the a priori restriction of the parameter space could play 
a crucial role in the identification of a particular isotope as critical. It is thus important 
to perform a general analysis in order to check whether the predictions obtained within the 
CQF for those nuclei close to a critical point are robust. We present below such an anal ysis 
in the region of the rare-earths. We follow closely the approach introduced in Ref. 
using catastrophe theory. In the next subsection the main ingredients of the theory are 
summarized and the relevant equations are particularized for the IBM Hamiltonian written 
in multipolar form, Eq. (^J). 
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A. The separatrix plane 



For the study of phase transitions in the IBM within the framework of catastrophe theory 
we aheady have the basic ingredients: the Hamiltonian of the system, Eq. ([Q), and the 
intrinsic state, Eq. (|12j). With them, we have generated the corresponding energy surface, 
Eq. ()14|) . in terms of the Hamiltonian parameters and the shape variables. It is our purpose 
to find the values of the parameters of the Hamiltonian that correspond to critical points. In 
principle this analysis involves the 6 parameters of the Hamiltonian, but a first simplification 
occurs since the energy surface only depends on 5 parameters : 

= 7^ + ff^M' + -^/^^-<3 7)+a3/3^ + |^), (15) 



where 



9 7 9 

e = Ed + 6 Hi - - + - + - 

4 5 5 
1 1 18 
ai = 4 ^0 + 2 '^^ 35 

02 = 2 V2 K,2 

uo = f . (16) 

Fortunately, it is possible to reduce the number of relevant (or essential) parameters to just 
two and stud y a ll phase transitions by using catastrophe theory jl^. We refer the reader 
to Refs. jl^-lisl for details of the application of this theory to the IBM case. The idea 
is to analyze the energy surface and obtain all equilibrium configurations, i.e. to find all 
the critical points of Eq. (|15p. First, the critical point of maximum degeneracy has to be 
identified. In our case, it corresponds to /? = 0. Next, the bifurcation and Maxwell sets 
are constructed Q|. Finally, the separatrix of the IBM is obtained by the union of 
Maxwell and bifurcation sets. In general a bifurcation set, corresponding to minima, limits 
an area where two minima in the energy surface coexist. A second order phase transition 
develops when these minima become the same. The crossing of a Maxwell set corresponding 
to minima leads to a first order phase transition. 

In order to follow this scheme, one has to identify the catastrophe germ of the IBM, which 
is the first term in the expansion of the energy surface around the critical point of maximum 
degeneracy that cannot be canceled by an arbitrary selection of parameters. In our case, 
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one finds that the first derivative in /5 = is always because of the critical character of the 
point for any value of the parameters. The second and third derivatives can also be canceled 
with an appropriate selection of parameters. However, if one imposes the cancellation of the 
fourth derivative, the energy becomes a constant for any value of (3. This means that the 
catastrophe germ is and the number of essential parameters is equal to two, which can 



be defined, following reference 



14 



3, 



as 



^ as-uo + e/jN -1) 
' 2a^ + e/{N-l)-as' ^ ' 

'^^ 2ai+£/(Ar-l)-a3' ^ ' 

where e, ai, 02, and 03 are defined in p6|) . The denominator in both expressions fixes the 
energy scale, which means that when it becomes negative, the energy surfaces are inverted. 
The essential parameters r\ and r2 can also be written in terms of the parameters appearing 
in ((H) as , 

g-(iV-l)(/.o-4/.2) , . 

1 £+(iV-l)(«o-3/€2 + i«4)' ^ ^ 

r = 4v^/.2(iV-l) 

£ + (iV-l)(/€o-3fi:2 + |fi:4)' ^ ^ 

A property of the parametrization used in this work is that the different chains of isotopes 
are located on a straight line that crosses the point corresponding to the f/(5) limit. The 
equation of this line is given by 

2ko — 7K2 + 11^4 

It should be remarked that the derivation of the essential parameters has nothing to do 
with catastrophe theory. The application of this theory begins once those parameters are 
obtained. The basic point is to translate every set of Hamiltonian parameters to the plane 
formed by the essential parameters ri and r2. This plane is divided into several sectors by 
the bifurcation set, that form the geometrical place in the parameter space where ^| = for 
a critical value of (3, and the Maxwell sets, the geometrical place in the space of parameters 
where two or more critical points are degenerate Q]. Both sets form the separatrix of the 
system, in this case of the IBM. In Ref. 14, Q| the IBM bifurcation (r2 axis, r2 = and 



ri < semi-axis, rn, and ri2) and Maxwell (negative ri semi- axis, rf^, and r-^^) sets were 
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obtained. They are all indicated in Fig. ^1 In this representation it is required that the 
denominator in (fT7j) and (fTHjl is positive. The separatrix for ri > is associated to minima 
while for ri < is associated to maxima (except the negative ri semi-axis). In order to 
clarify the figure on the separatrix, the energy surfaces corresponding to each set are plotted 
as insets. The half plane with r2 > corresponds to prolate nuclei, while the one with r2 < 
corresponds to oblate nuclei. Note that expressions ()19p and (^01) are only valid for prolate 
nuclei, but can be readily obtained for the oblate case. On this figure the symmetry limits 
and the correspondence with Casten's triangle are also represented. For completeness 
one should consider the case where the denominator of ()17|1 and (fTH|l is negative. It implies 
that the energy scale becomes negative and the energy surface should be inverted. The 
separatrix for this case is plotted in figure El and corresponds to the inversion of figure 
IT^ Again the schematic energy surfaces corresponding to each branch of the separatrix are 
shown as insets. Note that in this case the symmetry limits do not appear in the figure 
because they correspond to positive denominators for ri and r2. In our analysis only prolate 
nuclei are considered, because of that a new figure. Fig. El is included. In this figure, the 
right panel corresponds to positive denominators for ri and r2 while the left panel shows the 
case of negative denominator for ri and r2. In the following we will follow the convention 
presented in this figure. 

A set of parameters in the Hamiltonian corresponds to a point in the separatrix plane. 
The location of the point in that plane provides the required information on its transitional 
phase character. As mentioned above, it follows that points located on a separatrix line 
correspond to critical points. Note that the dynamical behavior of the system is controlled 
by the lowest minimum in the energy surface. In this sense we are adopting the Maxwell 
convention in the catastrophe theory language 3] and the only relevant branches of the 
separatrix are rjt^ and r2 = with ri < 0. All these branches correspond to first order 
phase transitions except for the single point (ri = 0, r2 = 0) that corresponds to a second 
order phase transition. The rest of Maxwell lines do not correspond to a phase transition 
because they are related to maxima. The interest of the bifurcation set, corresponding to 
minima, arises from the fact that it defines regions where two minima exits. In the following 
subsection the transitional isotope chains studied in this paper are analyzed in the separatrix 
plane. 
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B. Rare-earth region on the separatrix plane 



The fits presented in Sect. IIIII provide the parameter sets given in Tables HI and ITTl for the 
four isotope chains studied in this paper. In this section we plot the corresponding sequences 
of points representing the isotopes in each chain on the separatrix plane. As can be observed 
in the previous tables all the parameters for each chain are fixed except the value of Sd that 
changes along the chain. 

In figure HSl the positions of the different isotopes in the chains studied are plotted in the 
separatrix plane. The interpretation of these lines is given in Fig. As mentioned above, 
all isotopes in a chain lie on a straight line. The lighter ones are close to the ?7(5) point 
(spherical shapes) while as the number of neutrons is increased the corresponding points 
get increasingly away. For the heavier isotopes of Gd, and Dy the denominator of ri and r2 
becomes negative, which means that the left panel in Fig. has to be used. 

The main feature we find is that some nuclei are close to the Maxwell set r^: the closest 
are ^^^Nd (boson number = 8) and ^^°Sm (boson number A^ = 9) and not far away ^^^Gd 
(boson number A^ = 10). This can be complemented with the image of the energy surfaces 
plotted in Fig. The energy surface for ^^^Nd and ^^°Sm are rather flat around /3 = 0. 
For ^^^Gd the situation is not so clear. For Dy there is no isotope close to the critical point. 
According to our calculations, the transition from spherical to deformed occurs between 
A^ = 11 and A^ = 12. The isotope ^^^Dy is close to the Maxwell set but in the left panel. 
In this situation there should be two degenerate maxima. This can be observed in the 
corresponding energy surface (boson number A^ = 15) in Fig. ^2 The isotopes ^^°Nd 
(AT = 9) and ^^^Sm (A^ = 10) (also can be included in this situation ^^^Gd (A^ = 11) and 
i58]3y (^]S[ = X3)) are close to the bifurcation set r2 axis. Again inspection of Fig. ITT] shows 
that the energy surfaces for these isotopes has a minimum for f3 > and a maximum at 
(3 = 0. In figure El we show an amplification of the critical area. 

In conclusion, from this global analysis we find that ^^^Nd, ^^°Sm, and (less clearly) 
^^^Gd, are close to criticality. These isotopes are quite close but do not exactly coincide 



with previously proposed critical nuclei ^^°Nd and ^^^Sm 20|, |2J|, where the quite basic 
criterion was the closeness of their low-lying excitation spectra and transition intensities 
with the X{5) values. 
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C. Prediction of critical points within CQF 



The CQF uses a simplified Hamiltonian with only three parameters. For the description 
of transitional nuclei from the U (5) to the SU (3) limits the parameters are allowed to vary 
nucleus by nucleus. The representation of such calculations in the separatrix plane shows 
that all isotopes in a chain are basically on top of the straight line connecting the U{5) 
point, (ri,r2) = (1,0), and the SU{3) point, (ri,r2) = (— 4/3, 4v^/3). Note that this point 
corresponds strictly to the SU (3) Casimir operator. However, a more general CQF SU{3) 
Hamiltonian still lies very close to the latter point. In general, the same happens in the U{5) 
and 0(6) points. This means that within this framework the exploration of only a limited 
area in the separatrix plane is allowed. If all isotopes in an isotopic chain are forced to be 
located on the line connecting the U{5) and SU (3) points, it follows that one will more often 
find an isotope close to the (unique) critical point. In the calculations presented here we 
have seen that within the general formalism this is not always the case. For example, for 
Dy we did not find an isotope close to a critical point. 

In previous systematic studies in the rare-earth region using the CQF formalism, Ref. [3^ 
and |25|, the corresponding energy surfaces were not presented. We have constructed them 
from the parameters given in those references and the results obtained are consistent with 
those given in the present work. In particular, ^^^Nd and ^^'^Sm seem to be closest to a 
critical point. 



V. CONCLUSIONS 



In this paper we have analyzed chains of isotopes in the rare-earth region. In these 
chains nuclei evolve from spherical to deformed shapes. We have performed an analysis 
of the corresponding shape transitions to look for possible nuclei at or close to a critical 
point. We have used the more general one- and two- body IBM Hamiltonian and generated 
energy surfaces using the coherent state formalism. We have then used catastrophe theory 
to classify phase transitions and to decide if a nucleus is close to criticality. 

The approach used to fix the Hamiltonian parameters leads to a very good global agree- 
ment with the experimental data corresponding to excitation energies, B(E2)'s and val- 
ues. In particular, an excellent agreement with the measured values is obtained, which 
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is considered a key observable to locate phase transitional regions. The analysis presented 
here is consistent with previous CQF studies in the same region. As a result we find that 
-^^^Nd and ^^°Sm are the best candidates to be critical, but we should remark that ^^°Nd 
and ^^^Sm are not far away from it. 

A possible new way of defining critical nuclei is based on the "critical symmetries" E{5) 
or X(5) js, 0]. The properties associated with these solutions allow the identification of 
critical points by comparing the experimental data with characteristic energy and transition 
rate ratios. Thus, it may be possible to decide whether a nucleus is critical by analyzing 
its spectrum and decay properties. A trickier question is whether a fiat energy surface can 
be truly associated to a given nucleus with energy ratios close to X{5). A clear example is 
^^^Sm: in section QV Bl we have shown that according to our study the IBM energy surface 
of this nucleus is not so fiat as expected from previous analyses, i. e. in our work it does not 
correspond to a critical point as suggested earlier. However, if the spectrum and transition 
rates are analyzed (see figure ITOl and table ITTTl ). this nucleus reproduces reasonably well 
the main X{5) features. We note that in the general IBM framework there is no unique 
spectrum associated to a given potential energy surface, as implied by equations (|T7jl an 
flTH|) . Catastrophe theory constitutes a definite criterion regarding this issue, but does not 
provide a measurable signature in itself. 

It seems clear that further work is required to find experimentally identifiable features 
which signal criticality in an unequivocal way. 

VI. ACKNOWLEDGEMENTS 

We acknowledge insightful discussions with R.F. Casten, F. lachello, K. Heyde, J. Jolie, 
P. Van Isacker, and P. Von Brentano. This work was supported in part by the Spanish DGI- 
CYT under projects number FPA2000-1592-C03-02 and BFM2002-03315 and by CONACYT 
(Mexico). 



[1] D.J. Rowe, C. Bahri, and W. Wijesundera, Phys. Rev. Lett. 80, 4394, (1998). 
[2] F. lachello, N.V. Zamfir, and R.F. Casten, Phys. Rev. Lett. 81, 1191, (1998). 
[3] R.F. Casten, D. Kusnezov, and N.V. Zamfir, Phys. Rev. Lett. 82, 5000, (1999) 

15 



[4] J. Jolie, P. Cejnar, and J. Dobes, Phys. Rev. C 60, 061303, (1999). 

[5] F. lachello, Phys. Rev. Lett. 85, 3580, (2000). 

[6] F. lachello, Phys. Rev. Lett. 87, 052502, (2001). 

[7] F. lachello and A. Arima. "The interacting boson model". Cambridge University Press, 
Cambridge, (1987). 

[8] J. Jolie, R.F. Casten, P. von Brentano and V. Werner, Phys. Rev. Lett. 87, 162501, (2001). 

[9] O. Castanos, A. Frank and P. Federman, Phys. Lett. 88B, 203, (1979). 

[10] O. Castafios, P. Federman, A. Frank, and S. Pittel, Nucl. Phys. A 379, 61, (1982). 

[11] A. Frank, Phys. Rev. C 39, 652, (1989). 

[12] A. Gomez, O. Castafios, and A. Frank, Nucl. Phys. A 589, 267, (1995). A. Gomez, 

O. Castanos, A. Frank, C.E. Alonso, and J.M. Arias, Nucl. Phys. A 594, 483, (1995). 

[13] D.D. Warner and R.F. Casten, Phys. Rev. Lett. 48, 1385, (1982). 

[14] E. Lopez-Moreno and O. Castafios, Phys. Rev. C 54, 2374, (1996). 

[15] E. Lopez-Moreno, PhD thesis, Universidad Nacional Autonoma de Mexico, 1998, unpublished. 

[16] R. Gilmore. "Catastrophe theory for scientists and engineers". Wiley, New York, (1981). 

[17] N.V. Zamfir et al, Phys. Rev. C 60, 054312 (1999). 

[18] R.F. Casten and N.V. Zamfir, Phys. Rev. Lett. 85, 3584, (2000). 

[19] T. Klug, A. Dcwald, V. Werner, P. von Brentano and R.F. Casten, Phys. Lett. B 495, 55 
(2000). 

[20] R.F. Casten and N.V. Zamfir, Phys. Rev. Lett. 87, 052503, (2001). 

[21] A. Frank, C.E. Alonso, and J.M. Arias, Phys. Rev. C 65, 014301, (2001). 

[22] N.V. Zamfir et al, Phys. Rev. C 65, 044325, (2002). 

[23] B.G. Bizzeti and A.M. Bizzeti-Sona, Phys. Rev. C 66, 031301, (2002). 

[24] R. Kriicken et al, Phys. Rev. Lett. 88, 232501, (2002) 

[25] R. Fossion, C. De Coster, J.E. Garci'a-Ramos, T. Werner, and K. Heyde, Nucl. Phys. A 697, 
703, (2002). 

[26] A. A. Sonzogni, Nucl. Data Sheets 93, 599, (2001). 

[27] L.K. Peker and J.K. Tuli, Nucl. Data Sheets 82, 187, (1997). 

[28] M.R. Bhat, Nucl. Data Sheets 89, 797, (2000). 

[29] E. Dermateosian and J.K Tuh, Nucl. Data Sheets 75, 827, (1995). 

[30] Agda Artna-Cohen, Nucl. Data Sheets 79, 1, (1996). 

16 



[31] C.W. Reich and R.G. Helmer, Nucl. Data Sheets 85, 171, (1998). 

[32] R.G. Helmer, Nucl. Data Sheets 65, 65, (1992). 

[33] R.G. Helmer, Nucl. Data Sheets 77, 471, (1996). 

[34] C.W. Reich, Nucl. Data Sheets 78, 547, (1996). 

[35] R.G. Helmer and C.W. Reich, Nucl. Data Sheets 87, 317, (1999). 

[36] Bahaj Singh, Nucl. Data Sheets 93, 243, (2001). 

[37] E.N. Shurshikov and N.V. Timofeeva, Nucl. Data Sheets 67, 45, (1992). 

[38] G. Audi and A.H. Wapstra, Nucl. Phys. A 595, 409, (1995). 

[39] W.-T. Chou, N.V. Zamfir, and R.F. Casten, Phys. Rev. C 56, 829, (1997). 

[40] J.N. Ginocchio and M.W. Kirson, Nucl. Phys. A 350, 31, (1980). 

[41] A.E.L. Dieperink and O. Scholten, Nucl. Phys. A 346, 125, (1980). 

[42] A.E.L. Dieperink, O. Scholten, and F. lachello, Phys. Rev. Lett. 44, 1747, (1980). 

[43] P. Cejnar and J. Jolie, Phys. Rev. E 61, 6237, (2000). 

[44] J. Dukelsky, G.G. Dussel, R.P.J. Perazzo, S.L. Reich, and H.M. Sofia, Nucl. Phys. A 425, 93, 
(1984). 

[45] P. Van Isacker and J.Q. Chen, Phys. Rev. C 24 684, (1981). 

[46] H.E. Standley, "Introduction to phase transitions and critical phenomena" , Oxford University 
Press, Oxford (1971). 



17 



TABLE I: Values of £d in the Hamiltonian (in keV) for each isotopic chain as a function of the 
neutron number. 



Element 


Neutron Number 
84 86 88 90 92 94 96 98 100 


6oNd 


1686.3 1606.7 1645.4 1602.9 1536.1 1595.9 


62Sm 


1427.3 1393.5 1289.3 1210.8 1158.6 1192.5 1312.2 1452.0 


64 Gd 


1479.3 1508.7 1409.0 1300.4 1221.5 1174.4 1162.0 1176.5 


66 Dy 


1558.8 1607.6 1562.4 1503.9 1461.0 1427.7 1413.4 1409.2 1443.1 



TABLE IL Rest of the parameters in the Hamiltonian and in the E2 transition operator. 



Isotopes 


A (MeV) B (MeV) 


(keV) 


Ki (keV) K2(keV) 


Ks (keV) 


K4 (keV) 


eeff (e • b) 


X 


144-154^1 


16.75 


-0.51 


83.753 


-13.928 -17.151 


-101.27 


-187.57 


0.119 


-1.43 


146-160O 

62^™ 


18.05 


-0.46 


53.209 


-11.267 -14.674 


-31.769 


-131.24 


0.119 


-1.69 


148-162 pi 1 


22.55 


-0.76 


45.207 


-7.932 -13.129 


-35.224 


-156.24 


0.110 


-1.77 


150-166t~) 


25.06 


-0.80 


38.651 


-6.416 -13.638 


-59.165 


-163.05 


0.103 


-1.60 



TABLE IIL Relevant transition rates for ^^^Sm (in w.u.). 

Exp. X{5) This work CQF^'^) 



B{E2 




.0+ 


144 


144 


128 


144 


B{E2 


4^- 


-.2+ 


209 


228 


193 


216 


B{E2 


6^-- 


-4+ 


245 


285 


215 


242 


B{E2 


8^- 


-.6+ 


285 


327 


218 


248 


B{E2 : 


10+- 


-8+ 


) 320 


376 


210 


242 


B{E2 


0+- 


.2+ 


33 


91 


53 


57 


B{E2 


2t- 


-^4+ 


19 


52 


14 


20 


B{E2 


2+ - 

•^2 


-.2+ 


6 


13 


5 


11 


B{E2 


2t- 


.0+ 


1 


3 





0.1 


B{E2 


4- 


-.6+ 


4 


40 


7 


14 


B{E2 


4- 


-4+ 


5 


9 


2 


8 


B{E2 


4^- 


.2+ 


1 


13 


n ° 


0.1 






(a) 


Followin 


gRef. 


a. 



18 




o 



o 






Ijtheo 


□ 


3| exp 




3 J theo 



144 



146 



148 



150 



152 



154 



A 



FIG. 1: Excitation energies of Nd isotopes. 
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FIG. 2: Excitation energies of Sm isotopes. 
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FIG. 3: Excitation energies of Gd isotopes. 
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FIG. 4: Excitation energies of Dy isotopes. 
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FIG. 5: B{E2) transition rates for Nd isotopes. 
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FIG. 8: B{E2) transition rates for Dy isotopes. 
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FIG. 9: S2n values for Nd, Sm, Gd, and Dy isotopes. 
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FIG. 10: Spectrum of ""^^^Sm: (a) experimental, (b) X(5) symmetry, (c) this work, and (d) 
CQF Q. 
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FIG. 11: Energy surfaces for the different chain of isotopes. 
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FIG. 15: Representation of isotopes in the separatrix plane (with % < 0).The numbers on the 
isotopes correspond to the number of bosons. 
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